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Motivation 


• The detailed chemistry & physics models currently in 
vogue for simulations of shock structure in 
hypersonic flows were developed in the 1980s 

• Since that time, Navier-Stokes and Schrodinger 
equation solvers are ~10 6 times more powerful, 
but we are still using the same crude chemistry & 
physics models 

• For the last few years, NASA has been trying to 
remedy this situation 

Currently supported by Hypersonic-Entry Descent & Landing 
project in NASA's Space Technology Mission Directorate 
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Outline 


I. Introduction 

What are the shortcomings of the old models? 

II. Methods for computing reaction rate coefficients from first 
principles 

• Potential energy surfaces (PES) 

Map out the interatomic forces. 

• Quasiclassical trajectories (QCT) 

Simulate collisions using classical mechanics & Monte Carlo 
sampling. 

III. Analysis 

• 0-d Master Equations model 

Determine "phenomenological" rate coefficients. 

IV. Conclusions 
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Introduction 

The Apollo heatshield was designed with large safety 
margins, because hypersonic phenomena were not well 
understood 

Early observations and aerothermodynamics calculations 
clearly showed that thermodynamic and chemical 
equilibrium were not maintained in the bow shock layer 
during reentry 

Non-equilibrium radiation overshoot was observed, 
(predominantly atomic lines of nitrogen and oxygen) 

During the 1980s an empirical "2-temperature" 
chemistry model was developed 

Most of the parameters were obtained from older shock 
tube experiments 
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2-Temperature Chemistry Model 

Partitioning of internal energy modes 

? 

'rot ~ ' ^ 'electron 'vib ~ ' electronic 

Landau-Teller model used for post-shock internal 
energy relaxation 

dE vib (T v ) _ [E vib (T) - 

E vib C Ty )] 

dt Ty 

Internal energy relaxation decoupled from 
"chemistry" 

T ave = JS X Tib 1 ' 5 (usually T ave = (Tx TJX) 

^cUssi^r TJ — A •exp(-E fl /RT a J or (BT ave c ) •exp(-E A /RT a J 

Developed mostly by Chul Park at NASA Ames 
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Non-Equilibrium Overshoot at the Shock 


Immediately after the 
shock T » 10,000 K and 
T V ib ~ 300 K 

- ForT ROT =T prompt dissoc- 
iation & atomic excitation 

Radiation overshoot 

Closer to surface T = T V , B 
~ 10,000 K 

Thermal equilibrium 

Maximum T ~ 20,000 K for 
Trot= T and ~ 60,000 for 
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Our Objective 


• Develop a "physics-based" hypersonics model for 
high-speed Earth entry (> 10 km/s), validated 
against experimental data 

• Determine rate coefficients for dissociation, 
exchange and energy transfer reactions 

— N 2 dissociation for Earth entry 

- CO dissociation for Mars entry 

• Try to avoid preconceived assumptions about 
temperature 

— We have succeeded for internal energy modes, but still 
assume Maxwellian distribution for translation 
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Methods - PES 

• With all the advances in computer hardware and 
quantum chemistry software, achieving chemical 
accuracy (~4 kJ/mol) is still out of reach for 4-atom 
cases (e.g., N 2 +N 2 ) 

— We chose to employ a hybrid approach for the N 4 PES 

— Multi-Reference Configuration Interaction (MRCI) for 
geometries where the N-N bonds were greatly stretched 

— Single-reference method (CCSD(T)) for other geometries 

— Accurate empirical diatomic potential function for N 2 

• Energy computed for 1000s of geometries and fit to 
an analytic function for use in QCT calculations 

• For the N 3 PES only the MRCI method was used and 
combined with the empirical diatomic potential 
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Potential Energy Surface for N 2 + N 

Collisions 



Analytic representation of quantum chemistry energies for a 
grid of N 3 geometries 


Exchange reaction barrier 


color key 

blue - low energy 
red - high energy 
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Relative energies 
accurate to ±2 kcal/mol 


N 3 potential well 
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Methods - QCT 

Hamilton's equations of motion (q,, Pj are coordinates 

and momenta; H is total energy) 

<Uu_dH_ dpi__ _dH_ _ d VpES 

dt dpi dt dqi dqi 

Simulate collisions like N + N 2 or N 2 + N 2 

- Initial conditions correspond to N 2 (v,J) with random 
orientation, rotation and vibration phase, impact parameter 
(b) and relative collision velocity (V R ) 

- State-to-state cross section obtained from a large sample of 
trajectories 

rbmax 

Sr ( yirji> = I Pr (V R , b)bdb 

Jq 

- State-to-state rate coefficient obtained by integration of 
Boltzmann-weighted cross section over collision energy (E R ) 
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State-specific Rate Coefficients 


• N 2 X 1 ^ has 9390 rovibrational levels with the Leroy 
potential - we have obtained state-specific rate 
coefficients for all of them 

• For inelastic collisions rate coefficients exhibit relatively 
large values for |zW| >1 

• Inclusion of exchange reactions enhances this effect 

• Largest contributions to dissociation rate coefficient come 
from rovibrational levels near and above the dissociation 
energy 

• Quasi-bound levels contribute ~40% to thermal rate 
coefficient for T > 10,000 K 
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Rate Coeff cm3/s 


Thermal Rate Coefficients 


• Thermal dissociation rate coefficients agree well with 
experimental data (Appleton) 
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Analysis — 0-d Master Equations 



Master Equation simulate N 2 
relaxation-dissociation 

- All state-specific reactions for 
N 2 + N dissociation, exchange 
and energy transfer included 
(9390 ro-vibration levels) 

- Start with 97% N 2 + 3% N at 
300 K and impose and 
maintain a high value of T 

Follow the evolution of the 
N 2 rovibrational population 
from t = 0 to 10' 3 s 


6 8 10 
Energy [eV] 



Panesi et al., J. Chem. Phys. 138, 044312 (2013) 
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N 2 Relaxation Times (x) 

• For N 2 + N collisions Exchange reactions speed up relaxation 
(decrease x v and x R ) 

— x v = x R for T > 10,000 K 

• Dissociation decreases x v and x R for T > 20,000 K 

• Empirical Millikan-White x v >10x larger than present 
values 


with/without exchange with/without dissociation 



QSS? 

• Quasi-Steady State behavior is observed in Master 
Equation results 

- T v & T r reach plateaus at values less than T for T > 10,000 K 

- "Global" dissociation rate coefficient also plateaus 
For high T, most dissociation occurs before QSS is realized 

• x v ~ x R at high-T evident from plot of T v & T R vs t 
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Dissociation Rate Coefficient for N 2 + N 

• Global and thermal dissociation rate coefficients (k D and k D *, 
respectively) in good agreement with experimental data 
(Appleton) 


•Neglect of exchange 
reaction lowers k D by 
factor of 2 
•Quasi-bound levels 
make negligible 
contribution to k D 


io 4 /t 
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What the Master Equation Study tells us 


• For the non-equilibrium overshoot region T * T rot 

T > 20,000 K, x y ~ x R and T^|g — T R qj — Tj n ^ is a good 
approximation 

— Strong coupling between internal energy relaxation and 
dissociation 

— Maximum in T is much higher because rotational modes 
are not excited 

• Under these conditions, the dissociation rate 
coefficient is 1/3 lower than the thermal one 
- Exchange reaction lowers x and increases k D 
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Conclusions 


Accurate potential energy surfaces have been 
determined for N 3 and N 4 based on ab initio 
quantum mechanics calculations 

- Used in Quasi-Classical Trajectory calculations to obtain 
state-to-state rate coefficients for collisional processes 
important for hypersonic reentry conditions 

We have determined a full set of state-to-state rate 
coefficients for N 2 + N collisions 

- Excellent agreement with experimental data for thermal 
equilibrium conditions 

- Enables us to carry out a detailed study of non-equilibrium 
relaxation (Master Equation calculations) 

- A major assumption of the 2-temperature model (T R = T) is 
shown to be invalid for these conditions 
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Conclusions (2) 


• Similar calculations for N 2 + N 2 collisions have been 
carried out 

— Good agreement is found for dissociation rate coefficients 
at thermal equilibrium 

- Calculations of non-equilibrium effects are on-going 

• These state-to-state rate coefficients are being used 
in the development of coarse-grain reduced 
dimensionality models (stay for the next talk) that 
can be used in CFD calculations 

More detailed studies of the shock structure are needed 
to better characterize the non-equilibrium region 

Entry Systems & Technology Division 19 



Supporting material 
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Thermal Dissociation Rate Coeff icient 


• N 2 + N 2 dissociation 
rate coeff icient 
agrees well with 
shock-tube results of 
Byron and Appleton 

• No evidence for 
curvature in fit to 
Arrhenius expression 


Legend 
■ Byron 
A Appleton 
O Hanson 
★ QCT 

Dashed Park model 
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Park' s 2-Temperature Model for N 2 Dissociation 


Park’ s 2-T empirical model to 
account forT vib <T 

- Empirical 

- T = T s y T ( l s ) 

1 ave 1 x 1 vib 

- Originally Park used s = 0.5 

- Current model uses s = 0.7 

Model cannot be adjusted to agre 
with QCT data for full range of (T, 

"^"vib) 

- Rate coefficients fall off too much 
when T vjb « T 

Differences with QCT values indica 
dissociation is more strongly 
dependent on rotation than 
vibration 


Comparison of N 2 dissociation rate 
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of exchange reaction 


Activation energy for exchange is much less than 
for dissociation 

— 45 vs 225 kcal/mol 

Results in much larger changes in rotation- 
vibration energy than do inelastic collisions 




Exchange barrier on the N 3 PES 

Blue line - QCT 

Red line - quantum mech (Wong) 

Dashed line - QCT (Capitelli empirical PES) 
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Comments On The Development of Chemistry Models for Hypersonic Atmospheric 
Entry 

Richard Jaffe 

Aerothermodynamics Branch 
Entry Systems Technology Division 
NASA Ames Research Center 
Moffett Field, CA 94035 

Development of High-Fidelity Physics-Based Models to describe hypersonic flight 
through the atmospheres of Earth and Mars is underway at NASA Ames Research 
Center. The goal is to construct chemistry models of the collisional and radiative 
processes that occur in the bow shock and boundary layers of spacecraft during 
atmospheric entry that are free of empiricism. In this talk I will discuss our 
philosophy and describe some of our progress. Topics to be covered include 
thermochemistry, internal energy relaxation, collisional dissociation and radiative 
emission and absorption. For this work we start by solving the Schrodinger 
equation to obtain accurate interaction potentials and radiative properties. Then we 
invoke classical mechanics to compute state-specific heavy particle collision cross 
sections and reaction rate coefficients. Finally, phenomenological rate coefficients 
and relaxation times are determined from master equation solutions. 



Exchange reaction barrier 
Energy transfer: 

N 2 (Vj, Jj) + N 2 (v i ’,J i ’) - N 2 (v f ,J f ) + N 2 ( V f ,J/) 
N 2 (v jl J i ) + N ^ N 2 (v f ,J f ) + N 

Dissociation: 

N 2 (v i ,J i ) + — N + N + N 2 (v f \J f ’) 

N 2 (Vj, Jj) + N 2 (v j ’,J j ’) ^N + N + N + N 
N 2 (Vj, Jj) + N^N + N + N 



blue - low energy 
red - high energy 




N 3 Potential Energy Surface for angle 0= 115° 


